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Abstract 


This  report  presents  the  accomplishments  of  AFOSR  AASERT 
Grant  F49620-97- 1-0453.  The  focus  of  the  grant  was  the  use  of  the 
Application  Visualization  System  (AVS),  developed  by  Advanced  Vi- 
sual  Systems,  Inc.,  for  animating  supersonic  turbulent  flows  computed 
using  Large  Eddy  Simulation  (LES).  Three  flow  configurations  are  vi¬ 
sualized:  a  turbulent  flat  plate  boundary  layer,  an  8°  compression 
corner  and  a  25°  compression  corner.  The  animations  of  the  flat  plate 
supersonic  turbulent  boundary  layer  illustrates  the  developing  process 
of  the  large  eddies.  Large  turbulent  streaky  structures  are  clearly  ob¬ 
served  in  the  instantaneous  contour  plots  of  the  vorticity  modulus  in 
agreement  with  experiment.  The  qualitative  validation  of  generating 
the  inflow  condition  by  our  rescaling  and  reintroduction  method  is 
shown  by  comparing  the  flowfield  at  the  downstream  and  inflow  sta¬ 
tions.  Flowfield  animations  for  the  8°  compression  corner  display  the 
effect  of  the  shock  wave  on  the  flow  downstream  the  corner.  The  in¬ 
teraction  of  the  shock  wave  with  the  boundary  layer  and  the  incoming 
flow  is  visualized  through  comparing  the  temperature,  streamwise  ve¬ 
locity  and  streamwise  vorticity  at  the  cross  sections  before  and  after 
the  shock  wave.  The  flowfield  of  the  25°  compression  corner  over  a  pe¬ 
riod  of  time  are  visualized  by  a  sequence  of  animated  contours  of  the 
streamwise  velocity  and  the  temperature,  respectively.  The  instan¬ 
taneous  flow  separation  zone  is  indicated  by  the  negative  streamwise 
velocity.  The  significant  changes  in  the  flowfield  caused  by  the  strong 
shock  wave  formed  around  the  corner  are  found. 
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1  Introduction 


The  accurate  prediction  of  compressible  turbulent  flow  plays  a  significant 
role  in  engineering  design  including,  for  example,  propulsion  and  aerody¬ 
namics.  Computational  fluid  dynamics  (CFD)  is  widely  recognized  as  an 
effective  method  of  flow  analysis  for  high  speed  vehicles.  The  use  of  CFD 
can  greatly  reduce  the  length  of  the  design  cycle  and  the  cost.  Conventional 
Reynolds- averaged  Navier-Stokes  (RANS)  methods  yield  accurate  predic¬ 
tions  for  those  flow  regions  which  exhibit  little  or  no  separation.  However, 
they  fail  to  provide  accurate  predictions  of  critical  engineering  data,  e.g .,  the 
mean  and  rms  fluctuating  surface  pressure  and  surface  heat  transfer,  mean 
surface  skin  friction  and  locations  of  primary  and  secondary  separation  in 
those  flows  involving  strong  viscous-inviscid  interaction  [1].  Large  Eddy 
Simulation  (LES),  an  alternative  to  RANS,  has  been  demonstrated  to  be 
both  a  useful  research  tool  for  understanding  the  physics  of  turbulence  and 
also  a  predictive  method  for  flows  of  engineering  interest.  In  LES,  the  influ¬ 
ence  of  the  unresolved  scales  is  simulated  using  a  subgrid-scale  (SGS)  model 
or  the  inherent  dissipation  of  the  numerical  scheme,  denoted  by  Monotone 
Integrated  Large  Eddy  Simulation  (MILES)  by  Boris  et  al  [2]. 

Large  Eddy  Simulations  of  three  flow  configurations  (a  supersonic  flat  plate 
turbulent  boundary  layer,  an  8°  and  a  25°  compression  corner  at  Mach  3 
and  Res  =  2  x  104)  have  been  visualized  using  the  Application  Visualiza¬ 
tion  System  (AVS)  developed  by  Advanced  Visual  Systems,  Inc.  AVS  is  the 
leading  3D,  multi-platform,  data  visualization  environment  and  provides  an 
ideal  framework  for  the  delivery  and  implementation  of  visual  applications. 
The  details  of  the  grid  system  and  the  numerical  method  are  presented  in 
[5,  6,  7].  The  visualization  provides  insight  and  understanding  into  the  com¬ 
plex  physics  of  these  flows  which  are  characterized  by  very  large  databases. 
We  have  successfully  applied  AVS  to  animate  the  turbulent  flowfields  of  the 
three  flow  configurations.  The  turbulence  characteristics  (e.#.,  the  turbu¬ 
lent  streaky  structures,  the  development  and  movement  of  the  large  eddies 
and  their  interaction  with  the  shock  waves)  are  demonstrated.  Five  video 
animations  have  been  produced.  The  contents  of  the  videos  are  described 
in  the  following  section. 
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Results 


Five  video  animations  were  produced  (see  Section  3.1).  They  are  classified 
into  three  different  flow  configurations.  The  contents  of  the  videos  are  de¬ 
scribed  below.  Copies  of  the  videos  were  forwarded  to  AFOSR  as  they  were 
produced.  Additional  copies  can  be  obtained  from  Prof.  Doyle  Knight. 

2.1  Flat  Plate  Boundary  Layer 

A  Mach  3  supersonic  turbulent  flat  plate  boundary  layer  at  Res  =  2  x  104 
is  presented  in  video  no.  3.  The  streamwise,  transverse  and  spanwise  dimen¬ 
sions  are  Lx  =  14.8(5,  Ly  =  3.4(5  and  Lz  =  4.4(5,  respectively,  where  6  is  the 
incoming  boundary  layer  thickness. 

The  video  illustrates  the  accuracy  of  our  dynamic  turbulent  inflow  boundary 
condition  obtained  by  an  instantaneous  rescaling  of  the  boundary  layer  at  a 
downstream  location  to  the  inflow  boundary.  The  development  of  the  turbu¬ 
lent  inflow  boundary  condition  was  a  major  landmark  of  our  research.  Prior 
to  the  development  of  this  inflow  condition,  LES  of  a  supersonic  turbulent 
boundary  layer  was  computationally  expensive  because  the  computational 
domain  had  to  extend  from  upstream  of  the  leading  edge  of  a  flat  plate  (i.e., 
where  the  inflow  boundary  condition  was  known,  namely,  uniform  flow). 
Our  inflow  condition  allows  the  supersonic  turbulent  boundary  layer  to  be 
generated  with  a  finite  inflow  boundary  layer  thickness.  Our  inflow  condi¬ 
tion  extends  to  compressible  flows  the  method  of  Lund  et  al  [3].  At  each 
time  step,  the  flowfield  at  a  selected  downstream  station  ( x  =  11.8(5,  where 
x  is  measured  from  the  inflow  boundary  and  <5  is  the  inflow  boundary  layer 
thickness)  is  rescaled  and  reintroduced  at  the  inflow  boundary.  The  video 
displays  animated  plots  of  the  instantaneous  kinetic  energy,  temperature 
and  vorticity  modulus  at  z  =  2.26  with  closeups  at  the  inflow  and  outflow. 
The  video  shows  that  the  inflow  boundary  condition  does  not  create  any 
spurious  features.  Comparison  of  quantitative  turbulence  statistics  with  ex¬ 
perimental  data  confirms  the  accuracy  of  the  inflow  boundary  condition  [5]. 

The  video  also  illustrates  the  turbulent  structures  in  the  flat  plate  boundary 
layer.  A  particularly  striking  feature  is  the  strong  coupling  of  the  tempera¬ 
ture  and  velocity  field  evident  in  the  animation  of  the  combined  temperature 
contours  with  the  velocity  vectors  at  a  cross  section  of  x  =  11.8(5.  The  tur¬ 
bulence  production  in  the  viscous  sublayer  is  clearly  seen  and  agrees  with 
experimental  flow  visualization  [4].  The  entrainment  by  the  large  eddies 
generates  intense  changes  in  the  static  temperature  field.  A  cross  sectional 
combined  animation  of  streamwise  vorticity  with  the  velocity  vector  high¬ 
lights  the  regions  of  turbulence  production.  The  contours  of  the  vorticity 
modulus  at  y  =  0.02(5  and  y  —  0.1(5  clearly  exhibit  the  sublayer  streaks  in 
agreement  with  experiment. 


2.2  Supersonic  8°  Compression  Corner 

An  8°  supersonic  compression  corner  at  Mach  3  and  Res  =  2  x  104  is  pre¬ 
sented  in  video  nos.  1,  2  and  4.  The  computational  domain  is  -6.05  < 
Lx  <  65  where  x  is  measured  from  the  corner.  The  height  and  width  of  the 
computational  domain  are  Ly  =  3.4(5  and  Lz  =  4.4(5. 

Four  animations  are  presented  for  the  flowfield  at  j z  =  2.2(5.  The  animations 
show  the  instantaneous  kinetic  energy,  density,  temperature  and  streamwise 
vorticity.  The  shock  boundary  layer  interaction  significantly  affects  the  flow 
turbulence  as  indicated  in  the  instantaneous  contours  of  the  streamwise 
vorticity,  temperature  and  streamwise  velocity  at  two  different  cross  sections 
( x  =  ±3(5). 

The  movement  of  the  shock  wave  and  its  interaction  with  the  incoming  flow 
is  visualized  by  several  3-D  animations.  These  combine  an  instantaneous 
isosurface  of  static  pressure  (which  define  the  instantaneous  A— shock)  with 
contours  of  streamwise  velocity  (which  define  the  instantaneous  separation 
region).  Additionally,  the  pressure  isosurface  is  colored  with  the  instanta¬ 
neous  contours  of  static  temperature.  These  animations  illustrate  the  distor¬ 
tion  of  the  A— shock  by  the  incoming  large  eddies.  The  size  spanwise  ripples 
of  the  A—  shock  are  in  good  agreement  with  experiment  [5].  Although  the 
mean  flow  is  unseparated,  regions  of  local  instantaneous  separation  appear 
periodically. 

2.3  Supersonic  25°  Compression  Corner 

A  25°  compression  corner  at  Mach  3  and  Res  =  2  x  104  is  presented  in  video 
no.  5.  The  computational  domain  is  —20 5  <  x  <  8<5  where  x  is  measured 
from  the  corner.  The  height  and  width  of  the  computational  domain  are 
Ly  =  3.45  and  Ly  —  2.25. 

The  animations  are  composed  of  a  slice  within  the  boundary  layer  at  con¬ 
stant  distance  above  the  wall.  When  the  supersonic  turbulent  flow  inter¬ 
acts  with  the  shock  wave  generated  by  the  compression  corner,  significant 
changes  in  the  flowfield  appear.  Animations  of  the  contours  of  the  instan¬ 
taneous  streamwise  velocity  and  temperature  show  the  enhancement  of  tur¬ 
bulent  mixing  due  to  the  shock  boundary  layer  interaction.  The  mean  flow 
is  separated  at  the  corner.  The  instantaneous  flow  separation  zone  differs 
substantially  from  the  mean  separation  region  due  to  the  significant  motion 
of  the  A— shock.  The  spanwise  ripples  of  the  A— shock  are  in  good  agreement 
with  experiment  [6]. 
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